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Abstract 

In a confining, renormalisable, Dyson- Schwinger equation model of two-flavour QCD 
we explore the chemical-potential dependence of the dressed-quark propagator, 
which provides a means of determining the behaviour of the chiral and deconfine- 
ment order parameters, and low-energy pion observables. We find coincident, first 
order deconfinement and chiral symmetry restoration transitions at \i c = 375 MeV. 
/tt is insensitive to [i until fi ~ po = 0.7 ^ c when it begins to increase rapidly, 
is weakly dependent on (i, decreasing slowly with p and reaching a minimum 6% 
less than its ll = value at jj, = fiQ. In a two-flavour free-quark gas at ll = li c the 
baryon number density would be approximately 3/>o, where po = 0.16 fm~ 3 ; while 
in such a gas at the density is po- 
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1. Introduction. The production of a quark-gluon plasma and the quantita- 
tive exploration of its properties is a primary goal of experimental studies of 
high energy heavy ion collisions. Present efforts explore the domain of nonzero 
baryon number density [1]. The theoretical study of finite baryon density in 
QCD begins with the inclusion of a chemical potential, li, which modifies the 
fermion piece of the Euclidean action: ■ d + m — > 'y ■ d — 74/i + m. Whether 
conclusions reached via this application of equilibrium statistical field theory 
can actually be explored in heavy ion collisions is an unresolved question. 
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With the inclusion of \i the fermion determinant acquires an explicit imag- 
inary part, in addition to those terms associated with axial anomalies. The 
/i 7^ QCD action being complex entails that the study of finite density is 
significantly more difficult than that of finite temperature, T. As an illustra- 
tion of this, whereas numerical simulations of finite-T lattice-QCD have been 
quite extensive [2], requiring only asymmetric lattices, there is currently no 
numerical simulation algorithm for finite density lattice-QCD. Contemporary 
studies using the quenched approximation; i.e., eliminating the fermion deter- 
minant, encounter a forbidden region, which begins at /i — m n /2 [3]. Since 
m w — > in the chiral limit this is a serious limitation, preventing a reliable 
study of chiral symmetry restoration, for example. 

At T = = \i the Dyson- Schwinger equations [DSEs] have been employed ex- 
tensively in the study of confinement and dynamical chiral symmetry break- 
ing [4]; and the calculation of hadron observables [5]. As the quark-gluon 
plasma is characterised by deconfinement and chiral symmetry restoration, 
they also provide a continuum tool for the exploration of the onset and prop- 
erties of this phase of QCD. The DSEs are a system of coupled integral equa- 
tions whose solutions, the n-point Schwinger functions, are the fully-dressed 
Euclidean propagators and vertices for the theory. Once all the Schwinger 
functions are known then the theory is solved. To arrive at a tractable prob- 
lem one must truncate the system at a given level. Truncations that preserve 
the global symmetries of a field theory are easy to implement [6]. Preserving 
the gauge symmetry is more difficult but progress is being made [7]. 

An often-used approach is to focus on the DSE for the dressed-quark propaga- 
tor, S(p): the 2-point quark Schwinger function, whose kernel is constructed 
from the dressed-gluon propagator, D^ v (k), and the dressed-quark-gluon ver- 
tex, T^(k,p). Choosing Ansatze for D^ v {k) and T^(k,p) one obtains a single 
integral equation whose solution provides information about quark confine- 
ment and dynamical chiral symmetry breaking (DCSB). When the Ansatze 
are based on sound premises, such as the extensive body of DSE studies relat- 
ing to the gluon propagator [4,8] and quark-gluon vertex [9], and the results 
are insensitive to qualitatively equivalent variations of these Ansatze, then the 
conclusions of such a study can be judged robust. 

This approach has been used in a study of deconfinement and chiral symme- 
try restoration in 2-flavour QCD at finite-T [10]. Therein the quark DSE was 
solved using the one-parameter model dressed-gluon propagator of Ref. [11], 
which provided a good description of n and p-meson observables at T = = /i, 
and a continuum order parameter for deconfinement was introduced. That 
study established the existence of coincident, second-order deconfinement and 
chiral symmetry restoration transitions at T = T c pa 150 MeV with a critical 
exponent (3 = 0.33 ± 0.03, which is consistent with that of the N = 4 Heisen- 
berg magnet: (3 H = 0.38 ± 0.01. This has been argued [12] to characterise the 
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universality class containing 2-flavour QCD. Both the pion mass, m^, and the 
pion leptonic decay constant, /„-, were insensitive to T until T « 0.7 T c . How- 
ever, as T — > T c , the pion mass increased substantially, as thermal fluctuations 
overwhelmed quark-antiquark attraction in the pseudoscalar channel, until, at 
T = T c , — > and there was no bound state. These results confirm those of 
contemporary numerical simulations of finite-T lattice-QCD [2,13]; and make 
interesting the exploration of this model at finite chemical potential. 

2. DSE-model of two-flavour QCD. In a Euclidean space formulation, 
with {7^,7,,} = 2S^ V and 7^ = 7^, the renormalised dressed-quark propagator 
at fj, 7^ takes the form: 

S(p) = -i^-pa A {p) -il4^ P crc(P) + cr B (p) , (1) 



where p = (p, uj p = p^ + i/i), and satisfies the DSE: 



S^p)- 1 = ■ pA{p) + i lA u p C{p) + B{p) (2) 
= z£ij • p + Z 2 (i^ujp + m hm ) + . (3) 

Here mb m is the Lagrangian current-quark bare mass and the regularised self 
energy is 

E'(p) = 27 • P ^(p) + i74^p ^c(p) + E' B (p) ; (4) 

A 

*^(p) = / fo 2 ^(P " 9) I tr [^5(9)^(9,10] , (5) 

where: .F = A, B, C; V A = -(Zfij ■ p/\p\ 2 ), V B = Z u V c = -(Z^/up); 
and / = / d 4 q/(27r) A represents mnemonically a translationally invariant 
regularisation of the integral, with A the regularisation mass-scale. Although 
not explicitly indicated, the solutions, o>, are functions only of \p\ 2 and uj 2 . 

In renormalising we require that 

S(p, ^p)~ l \^ +p2=c2 =rf-p + ^4 + m R (() , (6) 



where ( is the renormalisation point and mn(() is the renormalised current- 
quark mass. This entails that the renormalisation constants are: 



^(c^A^) = l-s^( P > 4 )|g 2 +p; 
z 2 (c 2 ,a 2 ) = i-e' c (p> 4 ; 



1/1=0 



>+pi=C 2 



(7) 
(8) 
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m R (C 2 ) = Z 2 m hm + S / B (p,p4)|fp 1 2° +p 2 =C : 



(9) 



and yields the renormalised self energies: 



F(j>, Up) = & + S^(p, u p ) 



z'Ap,P4)\\ P r +P i=c- 



(10) 



where T — A, B, C; £a — 1 = £c and 



m R (( 2 ). 



In studying confinement the p-dependence of A and C is qualitatively impor- 
tant since it can conspire with that of B to eliminate free-particle poles in 
the dressed-quark propagator [14]. Furthermore, in the study of Ref. [15], the 
p-dependence of A and C was a crucial factor in determining the behaviour 
of bulk thermodynamic quantities such as the pressure and entropy; being 
responsible for these quantities reaching their respective Stefan-Boltzmann 
limits only for very large values of T and /i. It is therefore important in any 
DSE study to retain A(p) and C(p), and their dependence on p. 

The DSE- model is specified by the Landau-gauge choice [11] 



The first term in (12) is an integrable, infrared singularity that provides long- 
range effects associated with confinement [8]. The second term ensures that, 
neglecting logarithmic corrections, the propagator has the correct perturba- 
tive behaviour at large spacelike- A; 2 . Requiring that the large-distance effects 
associated with S 4 (k) are completely cancelled at small-distances by the sec- 
ond term fixes the ratio of the coefficients of these two terms in (12). Since 
Q(k 2 )/k 2 doesn't have a Lehmann representation it can be interpreted as de- 
scribing a confined gluon because the nonexistence of a Lehmann represen- 
tation is sufficient to ensure the absence of gluon production thresholds in 
iS-matrix elements describing colour- singlet to singlet transitions [4,17]. 

This model, (12), has no explicit yU-dependence, which can arise through quark 
vacuum polarisation insertions. As such it may be inadequate at large values 
of [/,, particularly near any critical chemical potential. However, until finite-/^ 
DSE studies of the gluon propagator become available no objective assessment 
of this possibility can be made. We therefore advocate proceeding with models 
such as (12) and assessing the qualitative and quantitative results obtained in 
the light of existing experiments and related theoretical studies. 





Q(k 2 ) 



4vr 2 d 47i 2 m 2 t S 4 (k) + 



1 _ e -[ fc2 /(4m t 2 )]- 



(12) 
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As in Ref. [11] we adopt the additional, simplifying specification: 

Fu(q,p)=7u, (13) 

which is often referred to as the "rainbow approximation". Using this trun- 
cation a mutually consistent constraint is: Z\ = Z 2 and = Z 2 [16]. In 
T = = /i studies this truncation is quantitatively reliable in Landau gauge; 
i.e., using (13) in (5) to obtain S(p) and calculating bound state masses via 
the Bethe-Salpeter equation, the choice of parameters in Q(k 2 ) that gives a 
good description is little modified by requiring the same quality fit with a 
more sophisticated vertex Ansatz, provided it is free of kinematic, light-cone 
singularities. Hence we do not expect the explicit form of the vertex to have 
a significant qualitative effect on our conclusions as long as it is in the class 
of Ansatze prescribed by Ref. [7]. This is supported by the results of Ref. [6], 
which indicate that a systematic improvement of (13) has very little effect on 
flavour-octet, pseudoscalar meson properties. 

In the present context, one can apply the confinement test introduced in 
Ref. [10] by analysing the configuration-space Schwinger function 

1 °° 

A 5 (r) = — / dp 4 e^ T a Bo (p = 0,u p ). (14) 



(Here and below the subscript "0" denotes a quantity calculated in the chiral 
limit, m R = [16].) 

As an illustration, consider a free, massive fermion, for which <Jb{p — 0, d p ) = 
Mj (ujp+M 2 ). This function has poles at p\ = — (M±//) 2 , which are associated 
with the //-induced offset of the particle and antiparticle zero-point energies, 
and one finds Ag(r) = 1/2 exp[— (M— /i) r] 6{M—y). This function is positive- 
definite and monotonically decreasing. Q 

In contrast, for a Schwinger function with complex-conjugate p 2 -poles, As(t) 
has zeros at r > [17]. Such a Schwinger function does not have a Lehmann 
representation and hence can be interpreted as describing a confined particle, 
with no associated asymptotic state. The confinement order parameter intro- 
duced in Ref. [10] is derived from this observation. Denoting by Tq the location 
of the first zero in A$(t), and defining k = 1/Vq, then deconfinement is ob- 
served if, for some \x = fi c , «(^c) = 0; at this point the baryon number density 



In the context of DSE studies the quantity M can be identified with a constituent- 
quark mass-scale characteristic of DCSB, which for n/ti-quarks is ~ 400 MeV. One 
therefore expects finite-// effects to become noticeable when /t ~ 400 MeV. 
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has overwhelmed the confinement mass-scale and the poles have migrated to 
the real axis. Obviously, in the case of a free particle k = 0. 



The simplest order parameter for DCSB is 

X = ^ 



So(|p| 2 = 0,a;J = V) 



(15) 



Its behaviour at any chiral symmetry restoration transition is identical to that 
of the chiral-limit vacuum quark condensate. 



We solve the DSE for S(p), (3), numerically with 

m t = 0.69GeV, m R {() = 1.1 MeV , 



(16) 



which were fixed in Ref. [11] by requiring a best x 2_ fit to a range of 7r-meson 
observables at T = = \i. We use d = 4/9 in (12) and renormalise at ( = 
9.47 GeV. With these choices our present study has no free parameters. 

The pion mass is given by [11] 



mX = W0(99) C ),; (17) 

A 

(m R (0(qq) ( )n = 8N c J B (a Bo - B [u 2 p a 2 c + \p\ 2 a\ + a 2 B ]) , (18) 
v 

which vanishes linearly with m R ((), and 



Nt = 2N r 



Bl{ 



p(T C (j' c + \p\ 2 (r A a' A + a B a' B 



(19) 



-! IpT (K 2 {?co»c ~ Wc?) + bT (°ao" a - W A f) + o B o" B - K) 2 ])} , 

with a' B = da b{\p\ 2 i / d\p\ 2 , etc. N w is the canonical normalisation constant 
for the Bethe-Salpeter amplitude in ladder approximation. The pion decay 
constant is obtained from 



U N n = 4N C J B [a A a B + § \pf (a' A a B - a A a' B )} 



(20) 



(17)-(20) were derived [11] under the assumption that the pion Bethe-Salpeter 
amplitude is given by T w = 275-EV The limitations of this assumption are dis- 
cussed in Ref. [16] but it is expected to be qualitatively reliable for small pion 
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Fig. 1. B(fi) from (21); B(p) > marks the domain of confinement and dynamical 
chiral symmetry breaking. The zero of B([i) is n c = 375 MeV. B(0) = (0.104 GeV) 4 . 

masses. The difference e = |1 — f n /N n \ is a measure of the error introduced by 
this representation; and e is no more than 0.17 over the range of fi considered. 
In solving the ladder Bethe-Salpeter equation with T w = 275-80, (17) provides 
an estimate of the pion mass obtained thereby that is accurate to within 1%. 

3. Results and Conclusions. It is now clear that in the domain of confine- 
ment and DCSB the solution of the DSE determines the order parameters for 
deconfmement and chiral symmetry restoration, and elementary pion observ- 
ables such as and /„.. Since T = jF(|p| 2 , u;^), they are all real. 

To explore the possibility of a phase transition we calculate the difference 
between the tree-level auxiliary-field effective-action [18] evaluated with the 
Wigner-Weyl solution, characterised by B = 0, and the Nambu-Goldstone 
solution, characterised by B ^ 0: 



B(jm) = 

A 



In 



(21) 



which defines a /i-dependent "bag constant" [19]. In (21), A and C repre- 
sent the solution of (3) obtained when Bq = 0; i.e., when dynamical chiral 
symmetry breaking is absent. This solution exists for all /i. 



B(n) is plotted in Fig. 1. It is positive when the Nambu-Goldstone phase is 
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Fig. 2. The order parameters for chiral symmetry restoration [x, diamonds] and 
deconfinement [re, circles]. fi c = 375 MeV. 

dynamically favoured; i.e., has the highest pressure, and becomes negative 
when the Wigner pressure becomes larger. The critical chemical potential is 
the zero of £>(/i); i.e., /i c = 375 MeV. This abrupt switch from the Nambu- 
Goldstone to the Wigner- Weyl phase signals a first order transition. The value 
of fi c in Ref. [15], obtained without the second term in (12); i.e., without the 
"perturbative tail" , is « 30% smaller. 

In Fig. 2 we plot the order parameters for chiral symmetry restoration, 
and deconfinement, re(/i), obtained from our DSE solutions. The chiral or- 
der parameter increases with increasing chemical potential up to /i c , with 
x(/ i c)/x(0) ~ 1-2, whereas re(/i) is insensitive to increasing /i. At fi c they both 
drop immediately and discontinuously to zero, as expected of a first-order 
phase transition. The increase of the chiral order parameter with fi is a nec- 
essary consequence of the momentum dependence of the scalar piece of the 
quark self energy, B(p), as is easily seen in Ref. [15] where the qualitative be- 
haviour of both quantities is identical. The vacuum quark condensate behaves 
in qualitatively the same manner as x- 

The behaviour of m n and f n is illustrated in Fig. 3. One observes that al- 
though the chiral order parameter increases with /z, m n decreases slowly as fi 
increases. This slow fall continues until fi 0.7 /x c , when m 7r ( y u)/m 7r (0) 0.94. 
At this point m w begins to increase although, for \i < /i c , m 7r (/x) does not 
exceed m 7r (0). This precludes pion condensation, in qualitative agreement 
with Ref. [20]. The behaviour of m n results from mutually compensating 
increases in ( 77^(0 (gg)^ and N%. This is a manifestation of the manner 
in which dynamical chiral symmetry breaking protects pseudoscalar meson 
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Fig. 3. Chemical potential dependence of the pion mass [m n , circles] and pion lep- 
tonic decay constant [/„., diamonds]. 

masses against rapid changes with /i. The pion leptonic decay constant is in- 
sensitive to the chemical potential until \x ~ 0.7 /i c , when it increases sharply 
so that fn((J>c)/ 'A (a* = 0) ~ 1.25. The relative insensitivity of and f w to 
changes in /i, until very near /z c , mirrors the behaviour of these observables 
at finite-T [10]. For example, it leads only to a 14% increase in the n — > /iv 
decay width at /i 0.9 fi c . The conjecture of Ref. [21] is inconsistent with the 
anticorrelation we observe between the /x-dependence of f n and m r 

We expect that improving upon the assumption = i^^Bq will only modify 
these observations to the extent that f n rises slowly and uniformly on < 
/i < /v The discussion of pion properties has relied implicitly upon the ladder 
Bethe-Salpeter equation. However, Ref. [6] indicates that improving upon this 
truncation; i.e., including additional skeleton diagrams in the quark DSE and 
pion Bethe-Salpeter equation in a manner that preserves Goldstone's theorem, 
will have little effect on pion properties and hence will not qualitatively affect 
our conclusions. 

The confined- quark vacuum consists of quark- ant iquark pairs correlated in a 
scalar condensate. Increasing /x increases the scalar density: (qq). However, as 
long as ji < n c there is no excess of particles over antiparticles in the vacuum 
and hence the baryon number density remains zero; i.e., V/x < /i c , p^~ d = [15]. 
This is just the statement that quark-antiquark pairs confined in the conden- 
sate do not contribute to the baryon number density. After deconfinement the 
quark pressure increases rapidly, as the condensate "breaks-up" , and an excess 
of quarks over antiquarks develops. At ji ~ 5/i c this quark pressure saturates 
the Stefan-Boltzmann limit: P u+ d = /i 4 /(27r 2 ). 
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As a gauge of the magnitude of /i c = 375 MeV we note that the baryon number 
density of a two-flavour free-quark gas at this chemical potential would be 

p u B F+dF (^) = ^ = 2.9 Po , (22) 
3 n z 

where po — 0.16 fm -3 is the equilibrium density of nuclear matter. For com- 
parison, the central core density expected in a 1.4 M & neutron star is 3.6- 
4.1 p [22]. Using this gauge, p corresponds to p ~ 260 MeV; i.e., p = 0.7 p c . 

Acknowledgments. We acknowledge useful conversations with D. Blaschke, 
A. H611, P. Maris and A. G. Williams. For their hospitality and support during 
visits where some of this work was conducted: AB and SS gratefully acknowl- 
edge the Physics Division at ANL; CDR, the Department of Physics at the 
University of Rostock, the Department of Physics and Mathematical Physics 
at the University of Adelaide during a term as a Distinguished Visiting Scholar 
and, in a subsequent visit, the Special Research Centre for the Subatomic 
Structure of Matter; and GIP, CDR and AWT, the Institute for Theoretical 
Nuclear Physics, University of Bonn, during the "Bonn Workshop on Confine- 
ment Physics" . This work was supported by: the Australian Research Council; 
the Department of Energy, Nuclear Physics Division, under contract no. W-31- 
109-ENG-38; Deutscher Akademischer Austauschdienst; the National Science 
Foundation under grant no. INT-9603385; and benefited from the resources of 
the National Energy Research Scientific Computing Center. 



References 



[1] T. D. Lee, Nucl. Phys. A 590 (1995) 11c. 

[2] J. Engels et al. Phys. Lett. B 396 (1997) 210. 

[3] M.-P Lombardo, J. B. Kogut and D. K. Sinclair, Phys. Rev. D 54 (1996) 2303. 

[4] C. D. Roberts and A. G. Williams, Prog. Part. Nucl. Phys. 33 (1994) 477. 

[5] P. C. Tandy, Prog. Part. Nucl. Phys. 39 (1997) 117. 

[6] A. Bender, C. D. Roberts and L. v. Smekal, Phys. Lett. B 380 (1996) 7. 

[7] A. Bashir, A. Kizilersu and M.R. Pennington, "The nonperturbative three point 
vertex in massless quenched QED and perturbation theory constraints" , e-print 
|hep-ph/9707421| . 



[8] See, for example, C. D. Roberts, "Confinement, Diquarks and Goldstone's 
theorem", e-print |nucl-th /9609039 , in Quark Confinement and the Hadron 
Spectrum, II (World Scientific, Singapore, 1997) 224; M. R. Pennington, 



10 



"Calculating hadronic properties in strong QCD", e-print hep-ph/9611242 ; and 
references therein. 

[9] See Ref. [7] and references therein. 



[10 

[11 
[12 

[13; 

[14 

[15; 

[16 

[17; 
[is; 

[19 
[20 
[21 

[22; 



A. Bender et al. Phys. Rev. Lett. 77 (1996) 3724. 

M. R. Frank and C. D. Roberts, Phys. Rev. C 53 (1996) 390. 

K. Rajagopal, "The Chiral Phase Transition in QCD: Critical Phenomena and 
Long- wavelength Pion Oscillations", e-print hep-ph/950431C . 

F. Karsch, Nucl. Phys. A 590 (1995) 367c. 

C. J. Burden, C. D. Roberts and A. G. Williams, Phys. Lett. B 285 (1992) 347. 

D. Blaschke, C. D. Roberts and S. Schmidt, "Thermodynamic properties of a 



simple, confining model" , e-print nucl-th/9706070 



P. Maris and C. D. Roberts, 'V- and X-meson Bethe-Salpeter amplitudes", 
e-print |nucl-th/9708029| , to appear in Phys. Rev C. 

P. Maris, Phys. Rev. D 52 6087 (1995). 

R. W. Haymaker, Riv. Nuovo Cim. 14 (1991) series 3, no. 8. 

R. T. Cahill and C. D. Roberts, Phys. Rev. D 32 (1985) 2419. 

H. Yabu, F. Myhrer and K. Kubodera, Phys. Rev. D 50 (1994) 3549. 

G. Brown, Nucl. Phys. A 488 (1988) 689c. 

R. B. Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C 38 (1988) 1010. 



11 



